p-Wave Resonant Bose Gas: A Finite-Momentum Spinor Superfluid 
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We show that a degenerate gas of two-species bosonic atoms interacting through a p- wave Feshbach 
resonance (as realized in, e.g., a 85 Rb- 87 Rb mixture) exhibits a finite-momentum atomic- molecular 
superfluid (AMSF), sandwiched by a molecular p-wave (orbital spinor) superfluid and by an s-wave 
atomic superfluid at large negative and positive detunings, respectively. The magnetic field can be 
used to tune the modulation wave vector of the AMSF state, as well as to drive quantum phase 
transitions in this rich system. 



A Feshbach resonance (FR) is an exceptionally fruit- 
ful experimental "knob" that allows exquisite tunability 
of interactions in degenerate atomic gases. This has led 
to realizations and studies of Bose-Einstein condensation 
(BECVBCS crossover of fermion-paired s-wave superflu- 
idity 0, 0, 0, 0|- The bosonic counterparts have also 
been extensively explored, and in fact, in the s-wave FR 
case, e.g., in 85 Rb [5J, predate recent fermionic develop- 
ments. As was recently emphasized 0, 0, Hj], in contrast 
to their fermionic analogs that undergo a smooth BEC- 
BCS crossover, resonant bosonic gases are predicted to 
exhibit magnetic field and/or temperature driven sharp 
phase transitions between distinct molecular and atomic 
superfluid phases. 

Motivated by these successes, recent attention has fo- 
cused on a realization of an even richer p-wave paired 
fermionic superfluidity 0, [l(J 11, [13], utilizing p-wave 
FR in K and 6 Li. Laboratory production of p-wave 
Feshbach molecules [ljl, [3] showed considerable promise 
toward this goal; however, reaching molecular degeneracy 
has been plagued with short molecular lifetimes. 

In another important development, experiments on a 
85 Rb- 87 Rb mixture have demonstrated a p-wave FR at 
B = 257. 8G between these two bosonic isotopes fljj . 
Although consequences of this two-body p-wave reso- 
nance on the degenerate state of such a gas mixture has 
not been further explored experimentally, it provides the 
main motivation for our work. In this Letter, we report 
on our study of a two-species degenerate Bose gas with 
a p-wave Feshbach resonant interspecies interaction. 

As summarized by the phase diagram in Fig. [TJ we 
find (within a mean-field treatment that we expect to 
largely survive fluctuations) that in addition to the nor- 
mal (N, i.e., nonsuperfluid) phase, the p-wave Feshbach 
resonant two-component balanced Bose gas (e.g., equal 
mixture of 85 Rb and 87 Rb atoms) exhibits three classes of 
superfluid phases: atomic (ASF), molecular (MSF), and 
atomic-molecular (AMSF) condensates. Our most inter- 
esting finding is that the AMSF, sandwiched between 
(large positive detuning) ASF and (large negative de- 
tuning) MSF phases is necessarily a finite-momentum Q 
spinor superfluid, with a characteristic wave vector (with 



h = i) 



< 



(i) 



tunable with a magnetic field (via FR detuning, v that 
primarily enters through the molecular condensate den- 
sity n m (v)), with a, m, £, and j p , respectively, the FR 
coupling, atomic mass, atom spacing, and a dimension- 
less measure of FR width|ll|. 




FIG. 1: (color online). Schematic temperature-detuning 
phase diagram for a two-species mixture of bosonic atoms 
exhibiting atomic (ASF), molecular (MSF), and atomic- 
molecular (AMSF) phases. In the AMSF p-wave, molecular 
condensate coexists with a finite-momentum Q (see left inset) 
atomic superfluid. The cartoon in AMSF phase illustrates the 
mechanism driving finite-momentum condensation. 

Within the narrow (jy p <C 1) FR approximation [|| 
0, [Hi 11 1, we find that the optimum collinear state is 
characterized by a single Q, Fulde-Ferrell-like state [l6j], 
as opposed to a +Q and — Q Larkin-Ovchinnikov- 
like sta tell7| found in unbalanced paired fermionic 



systems ljl [2(|. However, a more detailed study is 
necessary to ascertain the precise set of Q's at which 
AMSF condensation takes place. 

The physical picture behind the finite-momentum 
AMSF formation is illustrated in Fig. [TJ At intermedi- 
ate detuning where atomic gap closes within the MSF 
state, p-wave molecules decay via FR into a pair of 
atoms, which (due to the p-wave nature of the molecules) 
are necessarily created at finite and opposite momenta, 
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±k, and therefore at low temperature form a finite- 
momentum atomic condensate, AMSF. The energetic 
cost (~ k 2 /2m) of a finite momentum atomic condensa- 
tion is balanced by the lowering of the energy (~ aky/n m ) 
through FR hybridization between closed-channel p-wave 
molecule and open-channel pair of atoms that is only pos- 
sible at finite atomic momentum k, giving Q in Eq. ([I]). 

As we detail below, in addition, we find that each of 
the three superfluid classes (ASF, MSF, AMSF) in turn 
consist of distinct phases selected by detuning, tempera- 
ture, and background s-wave scattering lengths, and dis- 
tinguished by the nature (ferromagnetic or polar) of the 
p-wave molecular condensate and/or which combination 
of the two types of atoms is Bose condensed. We ex- 
plore the nature of these SF phases and associated phase 
transitions. 

To outline the derivation of these results, we consider 
a model of a gas mixture of two distinguishable bosonic 
atoms (e.g., 85 Rb, 87 Rb) [Hj], created by field opera- 
tors Vff(r) = (V'i(r), "02( r )) > an( i interacting through 
a p-wave Feshbach resonance associated with a tunable 
closed-channel bound state. The corresponding p-wave 
(I = 1) closed channel hetero-molecule (e.g., 85 Rb- 87 Rb) 
is created by a vector field operator ^(r) = (<%., </>t, 4>\). 
This system is governed by a grand-canonical Hamilto- 
nian H\tp a ,<f)\ = J d 3 rH, with 



(2) 



£7=1,2 



H.c. 



where 



with the 



molecular chemical potential p m = fi% +/12 — v adjustable 
by detuning v. We have taken atomic masses to be 
identical (a good approximation for the 85 Rb- 8 Rb mix- 
ture) [2l| . and focus on the balanced case of pi = ^2 = 
with [i fixing the total number of 85 Rb and 87 Rb atoms, 
whether in the (open-channel) atomic or (closed-channel) 
molecular form. 

The FR interaction encodes a coherent interconversion 
between a pair of open-channel atoms 1,2 (in a singlet 
combination of 1, 2 labels, as required by bosonic statis- 
tics) and a closed-channel p-wave molecule, with ampli- 
tude a 2l| . We have focused on a rotationally invariant 
FR interaction, with ui and a independent of the molecu- 
lar component i. This is an approximation to 85 Rb- 87 Rb 
mixture, where the p-wave FR around B — 257. 8G is split 
into a doublet by AB = 0.6G, similarly to the fermionic 
case of 40 k[I3, M, M, S3]. We leave the more realistic, 
richer case for future studies [13]. 

The FR coupling a and detuning v are fixed experi- 
mentally through measurements of the low-energy p-wave 
scattering amplitude[T3, [53] f p (k) = _ v -i + f kak 2^ lk - 3 - 
The scattering volume, v (tunable via v) and the charac- 
teristic wave vector, ko (a p-wave analog of the effective 



range, negative for the FR case) [ll| define our model 
parameters in terms of these experimental observables. 
The background (nonresonant) interaction density 

Ti-bg =7~la+ Ti-rn + 'Ham Consists of 

u - = E ^tm+^m^x, (3) 



|(0t.0)2 + ||^.^, 



(4) 
(5) 



a =1,2 



where coupling constants A CT , A12, 31,2, 9am are related 
to the corresponding s-wave scattering lengths (a±, 
etc.) in a standard way, and thus are fixed experimentally 
through measurements on the gas in a dilute limit. The 
miscibility of a two-component atomic gas requires 24 
a\a-i > a 2 2 , which may be problematic for the case of 
85 Rb- 87 Rb due to the negative background scattering 
length of 85 Rb. 

The molecular interaction couplings g±, gi (set by 
the L — and L = 2 channels of p-wave molecule- 
molecule scattering) and g am can be derived from a com- 
bination of s-wave atom-atom (A CT ) and p-wave FR (a) 
interactions [23j. 

Qualitative features of the phase diagram for the sys- 
tem can be mapped out through a mean-field treatment 
of the Hamiltonian This amounts to a minimization 
of the Landau free-energy functional F[^ a , <§] of classical 
fields ^(r), *(r)) corresponding to the coherent state 
field configurations for the atomic and molecular opera- 
tors. In the simplest approximation, F[^ a , <fr] takes the 
form identical to H[i/j a , cf)}, with the effective couplings 
(Mo-, /2 m , A (j , . . ■) that are functions of microscopic param- 
eters (p CT , v, A CT , . . .) appearing in 

Minimization of F\^ a , #] is quite straightforward 23] . 
For large positive detuning v, closed-channel molecules 
are gapped, and the ground state is a molecular vacuum. 
Thus, at low temperature p, m < and F^u,^ 
minimized by $ = 0, reducing to i^f*^] = F[^ ai 



fd 3 r [£ct=i,2 



+ Ai 2 |*i| 2 |* 2 | 



This functional is a special (Z7 (1) x U(l)) case of 
a O(N) x 0(M) model that has been studied 
extensively [53, [13, [13] ■ The free-energy is clearly 
minimized by uniform \E'i,\E , 2, as this lowers atomic 
kinetic energy. For A1A2 > A 2 2 in addition to the normal 
(nonsuperfluid) state, the system exhibits three ASF 
phases: (i) ASFi with *i ^ 0, * 2 = 0, (ii) ASF 2 with 
*i =0,4-2^ 0, (iii) ASF12 with ^ 0, * 2 ^ 0, sepa- 
rated by continuous phase transitions. For a balanced 
mixture pi = P2, the system exhibits a direct N-ASF12 
transition through a tetracritical point, pi = fa = 0, 
that is believed to be in the decoupled universality 
class [H El, S3]. For A1A2 < A 2 2 , the ASF i2 phase 
is absent, and ASFi and ASF2 are separated by a 
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that terminates at a bicritical 
All other transitions (N-ASFi, N- 



first-order transition 
point [H, [H, [2 
ASF 2 , and ASF r ASF 12 ) are in the XY universality 
class, breaking associated U(l) symmetries. The phase 
boundaries and the values of the atomic condensate 
order parameters can be straightforwardly computed 
within mean-field theory (MFT)[23], but are modified 
by fluctuations [lEHl- 

Within ASF phases, the spectrum of fluctuations can 
be straightforwardly computed by a Bogoliubov diago- 
nalization of coupled atomic and molecular excitations, 
with details depending on which of the three possible 
ASF phases is studied. In general, there will be one Bo- 
goliubov sound mode per broken atomic U(X) symmetry, 
with one Goldstone mode in ASFi and ASF 2 phases and 
two in ASF 12 [1|. 

In the opposite limit of large negative detuning, v open- 
channel atoms are gapped, and the ground state is an 
atomic vacuum. Hence, at low temperature ft < and 
F[^ ff , <&] is minimized by „ — and a uniform molecu- 
lar condensate <&, reducing to F m [<f>]/V = F[0, $\/V — 
-ft m \$\ 2 + f (** • *) 2 + f I* • $| 2 , for (orbital 1=1) 
spin-1 molecular bosons with well-studied thermodynam- 
ics H, 2^, 3(| 31 1 . In particular, we predict our system 
to also exhibit polar (MSF p for g 2 < 0) and ferromag- 
netic (MSFf m for g2 > 0) molecular condensates, respec- 
tively, corresponding to <fr = <I>on G [S2 x £/jv(1)]/^2 and 
<I> = $o(n + irn)j\pi S SO(3) order parameters, with 
n,rh,£ = nxrh an orthonormal triad and <!>o a complex 
amplitude, breaking £0(3) x CM 1 )© ill- The finite 
T N-MSF transitions are in the universality class of a 
complex 0(3) model(27j. 

The Goldstone mode content of these orbital I = 1 
molecular condensate phases is also identical to that 
of spinor condensates [28[, with MSF p exhibiting three 



E 



(MSF P 



(k) ~ k Bogoliubov modes and MSFf m char- 



acterized by one E, 



(MSF fn 



E 



(MSF fm ; 
ml 



(k) ~ k Bogoliubov and one 



(k) ~ k 2 ferromagnetic spin-wave modes. An 
attractive new feature of these orbital molecular con- 
densates is that g2 {y) is a tunable function of detun- 
ing that can be used to drive a 1st order MSF p -MSFf n 
transition 



11, 2 



To calculate the spectrum of low-energy excitations in- 
side the MSF (polar and ferromagnetic) phases, we sep- 
arate molecular field = + (p into a condensate and 
small fluctuations about it, obtaining H[ipa,^ + <p] ~ 



U)t\ where £f'[$] = W[0, #] is the 
zcroth-order approximation to the MSF ground state en- 
ergy, 



n 



(2) _ 



E 

(7 = 1,2 



■tl>le a ilj <7 +a$-i> 1 (-iV)ip2 + H.c. (6) 



is the quadratic atomic Hamiltonian density, and Ti 



(2) 



2 



€>*<£>• (p + H.c. is the 



quadratic molecular Hamiltonian density, with e„ = e a + 
g am \^\ 2 and w y = (u> + 5i|*| 2 )% + gi®*j<f>i + 2p 2 $*$j- 
To this quadratic order, the molecular and atomic ex- 
citations decouple and can therefore be diagonalized in- 
dependently. The molecular part has been extensi vely 
studied in the context of spinor F = 1 condensates |28|. 
For the polar (g 2 < 0) MSF p state, there are three 
"sound" modes, one Bogoliubov type with sound velocity 
c (msf p ) _ ^/((^ + g 2 )n m /2m, and other two degenerate 

spin- waves with velocity c^ 81 ^ = yj\g2\n m /2m. For the 
ferromagnetic {g% > 0) MSFf m state, there is one Bogoli- 
ubov mode, with sound velocity c ( MSFfm ) = \J g\n rn /2m, 
one quadratic k 2 /2m ferromagnetic spin- wave mode, and 
one gapped quadratic k 2 /2m + 2<7 2 n m canonically conju- 
gate mode. 

The atomic sector can also be readily diagonalized, giv- 
ing 



E, 



(MSF) 
(a) 



(k) = yfi 



* 2 I* -kl 2 , 



(7) 



where = k 2 /2m — /j, + g a mn m and n m is the molecular 
condensate density. The details of the spectrum only dif- 
fer quantitatively between the MSF p and MSFf m phases, 
both exhibiting a minimum at a finite k m j n for v > v*, 
and an atomic gap A(z/) = (& m i n ) that closes at 

the transition v Cl into the corresponding AMSF states. 
Such spectrum (displaying a finite k minimum) should be 
experimentally observable through, e.g., rf spectroscopy 
and is also a complementary way to detect the approach- 
ing finite Q instability (phase transition) into the AMSF 
state. Simple analysis inside MSF p gives 

, (MSFp) 



mm 
A (MSF P ) 

(MSFp) 

(MSFp-AMSFp) 



(2m 2 a 2 + mg)i 



mis, 



y^—(ma 2 + g)ma 2 n 2 n 
-(2ma 2 + g)n m , 
-(ma 2 + g)n m , 



ma^vn r , 



(8) 



with g = gi + g 2 — 2g am and we used lowest order MSF p 
relation fi m = 2fi — v w (g\ + g 2 )n m to eliminate the 
atomic fi in favor of molecular condensate n m and de- 
tuning v. The corresponding expressions inside MSFf m 
differ only slightly [23|], except that k^^ Ffm ^ lies in the n- 
rh plane perpendicular to the ferromagnetic quantization 
axis, t rather than along it as in the MSF p state. 

Upon further increase of v, the atomic gap A( MSF ) [y) 
closes at u ci and atoms Bose-condense at a fi- 

(MSF fm ) 
min 



k^ n SF) K) 



/ k (MSF p ) 
V min 



maJn mi kl 



nite 

ma^Jn m / V2), thereby breaking the remaining Z 2 x 
C^ajv(I) symmetry, with Z2 corresponding to the 
discrete part of (atom number) J7jv(1) unbroken 
in the paired MSF states. The two associ- 
ated order parameters are given by v P±(r) = 



E Q „ (±*Q»,ie 
the phase of Aq 



»(Q n T-e Q ) 4 
a* Q 



-Q, 



2 gjQ„ r\ w ]j ere Qq is 

• Q|e i9 c! an d Q n = n Q 
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(n G Z, an integer). The critical transition point for 
'$>+ and is split by ±|Aq|, respectively, so only 
condenses at v ci , allowing for the possibility of two tran- 
sitions, MSF-AMSF_ followed by AMSF_-AMSF + . The 
values of order parameters and transition points can be 
straightforwardly worked out. Within MFT^we find that 
a single-Q state is energetically preferred [23J, but we do 
not expect this to survive generically. 

In addition to the molecular and atomic superfluid- 
ity, AMSFP' _ phases arc finite-momentum condensates 
that generically can exhibit crystalline order and there- 
fore are supersolids (at least in the case of more than one 
Q„ condensation [34j ) . For a collinear set of Q n 's, the su- 
persolid is a unidirectional density wave breaking trans- 
lational invariance along Q n , with the latter aligned with 
the quantization axis, n of the MSF p state and transverse 
to the quantization axis I (i.e., lying in the fi-rh plane) of 
the MSFf m state. Thus, while rotational 0(2) symmetry 
about the n axis remains intact inside the AMSF p state, 
it is spontaneously broken by such a uniaxial density wave 
inside the AMSFf m state. It is notable that in this lat- 
ter case, in the presence of fluctuations, such superfluid 
density wave will exhibit quantum liquid-crystal phe- 
nomenology similar to that of the fermion-paired Larkin- 
Ovchinnikov superfluidfH, HI]. The excitation spectra 
inside AMSF phases can be computed via a generalized 
Bogoliubov transformation and involve a diagonalization 
of a 10 x 10 matrix (corresponding to three and one cou- 
pled complex molecular and atomic fields), leading to 
Goldstone modes consistent with above symmetry-based 
arguments [2^1 . 

The nature of MSF Pi f m -AMSF P; f m transitions (beyond 
MFT) remains an open question. Based on the experi- 
ence with spinor condensates [32] and LO superfluid [3^], 
we expect this system to exhibit a variety of fractional 
composite topological defects. We leave detailed study of 
these to a future publication^. 

Upon further increase of v, the AMSF transitions into 
the ASF at v c2 sa (2\-g am )n a (with A = ±(A 1 +A 2 +2A 12 ) 
and the atomic condensate density n ), determined by 
the point of vanishing of the n m {v) and Q{v). The AMSF 
can thus exhibit a broad stability range z/ c2 — v c \ , set by 
a combination of n a , n m , \i, gi, that in 3d we expect 
to survive beyond our mean- field analysis. 

To summarize, we studied a degenerate gas of two- 
species bosonic atoms interacting through a p-wave Fes- 
hbach resonance, as realized in a 85 Rb- 87 Rb mixture. We 
showed that at intermediate FR detuning, such gas ex- 
hibits an atomic-molecular superfluid (AMSF) state con- 
densed at a finite momentum, that undergoes phase tran- 
sitions into a molecular p- wave (orbital spinor) superfluid 
(MSF) and into an s-wave atomic superfluid (ASF) at 
large negative and positive detunings, respectively. A 
magnetic field can be used to tune the modulation wave 
vector of the AMSF between zero and a value set by in- 
teractions as well as to drive quantum phase transitions 



in this rich system. 
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